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The differential rotation of e Eri from MOST data* 
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From high-precision MOST photometry spanning 35 days the existence of two spots rotating with slightly differing periods 
is confirmed. From the marginal probability distribution of the derived differential rotation parameter k its expectation 
value as well as confidence limits are computed directly from the data. The result depends on the assumed range in 
inclination i, not on the shape of the prior distributions. Two cases have been considered: (a) The priors for angles, 
inclination i of the star and spot latitudes /3i,2, are assumed to be constant over i, /3i, and /32. (b) The priors are assumed 
to be constant over cos i, sin f3i, and sin 02. In both cases the full range of inclination is considered: 0° < i < 90°. Scale- 
free parameters, i. e. periods and spot areas (in case of small spots) are taken logarithmically. Irrespective of the shape of 
the prior, k is restricted to 0.03 < k < 0.10 (one-<r limits). The inclination i of the star is photometrically ill-defined. 
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1 Introduction 

In late-type stars with their deep convection zones non-uni- 
form rotation is driven by the action of the Coriolis force 
on that convective turbulence (cf. Kitchatinov & Riidiger 
1993). Now quantitative models of differential rotation for 
the Sun and solar-like stars are available (Riidiger & Kitcha- 
tinov 2005; Kiiker & Riidiger 2005; Kiiker [these proceed- 
ings]) which should be compared with real stars. 

The outcome of these theoretical efforts is a lapping 
time which is for a given main-sequence star nearly inde- 
pendent of its rotation period (in the case of the Sun, a G2 
dwarf, roughly 100 days). 

Precision photometry of a spotted star with spots differ- 
ing in latitude allows a direct measurement of the differen- 
tial rotation parameter k. It parameterizes the surface rota- 
tion. With P cq denoting the equatorial rotation period that 
at latitude (3 is Pp = P cq /(l - k sin 2 0). 

Here the MOST data (Croll et al. 2006; Croll 2006) of 
the star e Eri have been reanalyzed in a Bayesian frame- 
work. The motivation was to get realistic error estimates by 
computing k's marginal distribution and to find out, how 
it depends on the chosen prior. In the essence a Bayesian 
approach explores the whole likelihood mountain in the N- 
dimensional parameter space. It looks not only to the most 
probable set of parameter values, but provides expectation 
values. Moreover, from a marginal distribution reliable con- 
fidence limits can be easily derived. 



* Based on data from the MOST satellite, a Canadian Space Agency 
mission, jointly operated by Dynacon Inc., the University of Toronto In- 
stitute for Aerospace Studies and the University of British Columbia, with 
the assistance of the University of Vienna. 
** Corresponding author: HEFroehlich@aip.de 
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Fig. 1 The e Eri light curve with an obvious trend removed. 



The nearby K2 dwarf is young (< 1 Gyr), rotates twice 
as fast as the Sun, and shows strong chromospheric activity 
(cf. Biazzo et al. 2007). 



2 The MOST data 



The MOST photometric satellite (Walker et al. 2003) ob- 
served three consecutive rotations of e Eri in 2005. The data 
consist of 492 data points, one point per orbit of the satellite. 
The point-to-point precision of the data is 50 ppm rms (i.e. 
± 0.00005 mag). The rectified light curve (Fig. 1) spanning 
35 days shows an overall variation of a hundredth of a mag- 
nitude. 
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3 A Bayesian data analysis 

In the case of two circular spots the light curve is deter- 
mined by nine parameters: two periods, two epochs, two 
latitudes, two areas, and the inclination of the star. The pa- 
rameters specifying limb darkening and spot rest intensity 
are assumed to be given (cf. Croll et al. 2006). 

There are nuisance parameters: an offset in the photo- 
metric zero point and a long term trend in the data. With 
these two parameters introduced shifting the light curve ver- 
tically and even adding a trend does not alter the results. The 
flux error is considered a Gaussian with unknown variance. 
Integrating away that error a, assuming Jeffreys 1/er prior, 
sounds strange, but perhaps it is not so bad an assumption 
and it leads to quite a simple formula for the likelihood. The 
likelihood function used is therefore a mean likelihood with 
respect to these three parameters. One can say the method 
itself determines offset, slope and measurement error. 

The result, especially the star's inclination i, depends 
on the prior distribution functions. In the following it is as- 
sumed that rotational frequencies as well as spot areas (if 
small compared with the star), i. e. parameters missing a 
characteristic scale, have constant priors if taken logarith- 
mically. Two cases are considered: In Case A the prior is 
assumed constant over the inclination i and the spot lati- 
tudes ft, and ft as in the Croll (2006) paper. In Case B it is 
assumed constant over cos i, sin ft, and sin ft, respectively. 
An inclination prior constant over cos i means that nothing 
is a priori known about the orientation of the rotational axis. 
A latitude prior constant over sin f3 takes into account that 
there are more possibilities to locate by chance a spot near 
the equator than near the poles. 

In order to explore the likelihood mountain over a nine- 
dimensional parameter space the Markov chain Monte Car- 
lo (MCMC) method (cf. Press et al., 2007) has been applied. 
The results should be therefore best compared with the anal- 
ysis of the MOST data by Croll (2006), who has employed 
that MCMC technique. 

A set of 64 Markov chains was generated. Each chain 
has performed 10 7 steps. After a burn-in period of 1000 
steps every 1000th successful step has been recorded to sup- 
press the correlation between successive steps. For model- 
ing light curves Budding's star-spot model (1977) has been 
used. 

4 Results 

Expectation values and modes as well as one-cr confidence 
limits are presented in Table 1, augmented by the solution 
of Croll (2006; his "Wide Prior" Case). The Case-B prob- 
ability distributions for the parameters k and i are shown 
in Figs. 2 and 3, respectively. Especially the inclination pa- 
rameter i proves ill-defined by photometry alone. There is 
no strong correlation between i and k (Fig. 4). 

The differences between the two cases A and B are mar- 
ginal, i. e. the outcome is rather insensitive to the shape of 
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Fig. 2 The marginal distribution for the differential rotation pa- 
rameter k (Case B). Vertical lines denote expectation value and 
the 68 per cent confidence region. Additionally to the marginal- 
ized distribution the run of mean likelihood (dashed) is given, pro- 
viding an impression how the goodness of fit varies with k. The 
mismatch hints at the non-Gaussianity of the likelihood mountain 
(cf. Lewis & Bridle 2002). 
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Fig. 3 The marginal distribution for the inclination i (Case B). 
Vertical lines denote expectation value and the 68 per cent con- 
fidence region. The high-i solution is suppressed by the chosen 
prior. 

the prior. Most important is the restriction of the inclination 
i to values between 15° and 40° in the "Wide Prior" Case 
of the Croll MCMC analysis. (The marginalized likelihood 
of i, the last plot of his Fig. 3, indicates that inclinations 
beyond 40° are not ruled out.) Here the whole range, 0° < 
i < 90°, is taken into account. It includes a second peak in 
the likelihood mountain: i w 72°, ft w 61°, ft « 73°. This 
high-i solution is even more probable, by a factor of 4.6, 
then the best low-i solution (i w 24°, ft « 14°, ft w 25°). 

Contrary to the Croll solution the second spot is always 
visible. 

By fitting spectroscopic measurements (2000/2001) to 
a spot/plage model Biazzo et al. (2007) find two spots but 
larger than ours. By the way, this could be due to the lower 
temperature contrast between spot and photosphere these 
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Table 1 Results: Given are the expectation values, the modes (Case B) and la errors. Periods P are in days, epochs are with regard to 
E = HJD - 2,451,545. Note that the Case- A i and f3\ are both bi-modal. To give only one interval, for once the 90% confidence region 
is chosen. 
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Fig. 4 The parameters of differential rotation k and inclination i 
are somewhat correlated (Case B). 



authors have found. Perhaps the spots are persistent. The 
stated spot latitudes, /?i w 21° and [3 2 ~ 48°, are within the 
uncertainties of our estimate. 

Equatorial rotational period (P oq = 11.2 d) and radius 
(0.72i? Q ) of e Eri are well known. So in principle spec- 
troscopic determinations of the projected rotational velocity 
may restrict the inclination i. Unfortunately, e Eri rotates 
slowly. Saar & Osten (1997) find v sini w 1.7 ± 0.3 km/s. 
This leads to a sin i = 0.5±0.3 hinting at the low-i solution, 
but is too uncertain as to constrain the range of i values very 
much. A recent compilation by Valenti & Fischer (2005) 
gives v sini w 2.4 ± 0.5 km/s. This favours an inclination 
close to 50°. 

One should note that the planetary companion (Hatzes 
et al. 2000) as well as a 130 AU dust ring show both an 
inclination distinctly below 30°, namely 26°. 2 (Benedict et 
al. 2006) and w 25° (Greaves et al. 2005), respectively. 



5 Conclusions 

As expected, with 0.03 < k < 0.10, the estimated differen- 
tial rotation parameter k proves smaller than that of the Sun 
(fc Q w 0.2). The horizontal shear is 0.017 < 5£l < 0.056 
rad/d, the lapping time 130 d. An independent estimate of 
the star's inclination would exclude either the low-i solution 
or the high-i one and would help to constrain k even more. 

Acknowledgments: The author thanks the MOST team 
for the e Eri data and W.W. Weiss from Vienna for valuable 
discussions. 
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